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We investigate the one-body dissipation dynamics in heavy-ion collisions of '*O+!°O using a fully 
three-dimensional time-dependent Hartree-Fock (TDHF) theory with the modern Skyrme energy 
functional and without any symmetry restrictions. The energy dissipation is revealed to decrease 
in deep-inelastic collisions of the light systems as the bombarding energy increases owing to the 
competition between collective motion and single-particle degrees of freedom. The role of spin-orbit 
force is given particular emphasis in deep-inelastic collisions. The spin-orbit force causes a significant 
enhancement of the dissipation. The time-even coupling of spin-orbit force plays a dominant role at 
low energies, while the influence of time-odd terms is notable at high energies. About 40-65% of the 
total dissipation depending on the different parameter sets is predicted to arise from the spin-orbit 
force. The theoretical fusion cross section has a reasonably good agreement with the experimental 
data, considering that no free parameters are adjusted to reaction dynamics in the TDHF approach. 


PACS numbers: 24.10.-i, 25.70.-z, 21.60.Jz 


I. INTRODUCTION 


Time-dependent Hartree-Fock (TDHF) theory, orig- 
inally proposed by Dirac [1], has been widely ap- 
plied to describe the fusion excitation functions, deep- 
inelastic collisions, fission, collective excitation, and nu- 
clear molecular resonances; for reviews see Refs. [2, 3]. 
It provides, from a fully microscopic point of view, the 
dynamical foundation of large amplitude collective mo- 
tion [4]. TDHF theory has exactly treated the one- 
body dissipation caused by the collisions of nucleons with 
mean-field potential and nucleon transfer between the 
colliding partners. The two-body dissipation, which was 
included in quantum molecular dynamics models [5-7] 
and is likely to play a significant role at high energies, 
has been neglected within the framework of mean-field 
dynamics. 

A number of approaches based on TDHF have been 
developed to investigate the dissipation mechanism in 
heavy-ion collisions and resonance dynamics. Density- 
constrained TDHF (DC-TDHF) [8] has utilized the time- 
dependent instantaneous densities as the constraints to 
perform the static Hartree-Fock minimization, thus ca- 
pable of extracting the excitation energy [9] and the un- 
derlying nucleus-nucleus interaction potential [8]. The 
method is the dynamical analog of static adiabatic ap- 
proximation and has been applied to calculate the fusion 
cross section at both sub- and above-barrier energies [10- 
12]. In the dissipative dynamics TDHF (DD-TDHF) ap- 
proach [13, 14], mean-field evolution was assumed to be 
properly reduced to one-dimensional dissipative dynam- 
ics. By using the macroscopic reduction procedure, the 
one-body dissipation in the entrance channel of heavy- 
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ion fusion reactions has been extracted from the simu- 
lations of microscopic TDHF [13, 14]. The technique of 
Wigner distribution function [15, 16] has been applied 
to probe the single-particle dissipation in TDHF from 
a phase-space perspective. As indicated by Loebl and 
collaborators, the one-body dissipation is not enough to 
achieve the true equilibrium as the reaction proceeds to 
longer contact times [15]. The tensor force has been in- 
cluded in recent TDHF calculations [17-19]. The time- 
even contribution of the tensor force was found to have 
important roles in heavy-ion collisions [17, 18]. The full 
tensor expression was introduced to study the resonance 
dynamics [19]. 


Since TDHF theory is based on the independent- 
particle approximation, the expectation value of one- 
body observables can be well described. However, the de- 
scription of dynamical fluctuation and two-body dissipa- 
tion requires the theoretical development beyond TDHF. 
The stochastic extension of TDHF theory [20, 21] has 
been proposed to include the dynamical fluctuations of 
collective variables beyond mean field in a fully micro- 
scopic framework. The Balian-Vénéroni variational prin- 
ciple [22, 23] was first applied to realistic calculations 
of fragment mass and charge distributions in heavy-ion 
collisions by Simenel [24]. By comparing with the stan- 
dard TDHF calculations, an increase of the fluctuations 
has been found and the agreement with experimental 
data has been much improved. The two-body dissipation 
arising from nucleon-nucleon collisions has been taken 
into account in a quantum mechanical way with the ex- 
tended TDHF [25, 26] or time-dependent density matrix 
(TDDM) approaches [27-29]. So far only a few studies 
on the dynamical fluctuation and two-body dissipation 
have been implemented because beyond TDHF calcula- 
tions require a lot of numerical efforts and computational 
time. 
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The dissipation dynamics has an influence on the be- 
havior of nucleus-nucleus collision. The energy dissipa- 
tion, closely associated with the interplay between the 
collective motion and the single-particle motion, depends 
on the delicate balance between reaction time and rear- 
rangement time of the mean field. When the collective 
motion is slow enough so that the mean field has enough 
time to rearrange itself, the reaction dynamics follows 
the adiabatic motion. For example, the experimentally 
observed Landau-Zener effect [30, 31] in inelastic cross 
section is due to the breakdown of adiabatic condition. 
The Landau-Zener effect and its applications to heavy- 
ion collisions were, for the first time, investigated with 
TDHF theory by constructing the boost-invariant adia- 
batic single-particle states [32]. 


Earlier TDHF calculations systematically underesti- 
mated the dissipation of energy from the relative col- 
lective motion into internal degrees of freedom due to 
the simplified effective interaction and symmetry restric- 
tions [33-36]. For example, there has been the puzzle of 
the fusion window anomaly, which was later solved by in- 
cluding the time-even terms of spin-orbit force in TDHF 
calculations [37]. The fusion excitation function was en- 
hanced by as much as 20% for the 4°Ca+*°Ca system by 
removing the isospin symmetry restriction on the nuclear 
wave functions [33]. In Ref. [38] it was demonstrated 
that, despite using various approximations, the time-even 
spin-orbit force has a significant effect on the dissipation 
dynamics in heavy-ion scattering. However, the time- 
even spin-orbit interaction and the three-dimensional ge- 
ometry have not been incorporated simultaneously in 
these earlier TDHF calculations, and also the time-odd 
terms of spin-orbit force have been neglected. The first 
three-dimensional TDHF calculations including the full 
spin-orbit interaction [39, 40] have been performed using 
the TDHF3D code. The role of spin-orbit interaction on 
single-particle states in heavy-ion collisions was discussed 
in Ref. [41]. The time-odd terms of spin-orbit force was 
addressed to be the origin of the spin excitation mecha- 
nism [42]. 


With the development of computer ability, a fully 
three-dimensional TDHF calculation with modern en- 
ergy functional and without any symmetry restrictions 
has been realized in recent years. It is expected to pro- 
vide a better description of heavy-ion collisions [43-50] 
and resonance dynamics [51-56]. The purpose of present 
article is to systematically investigate the dissipation dy- 
namics in deep-inelastic collisions and the role of spin- 
orbit force by using modern TDHF calculations. The ar- 
ticle is organized as follows. In Sec. II, we briefly explain 
TDHF theory and the details in numerical calculations. 
In Sec. III, we illustrate the dissipation dynamics and 
the role of time-even and time-odd terms of spin-orbit 
force. The fusion cross section obtained from modern 
TDHF calculations is compared with experimental data. 
A summary is given in Sec. IV. 


I. THE TIME-DEPENDENT HARTREE-FOCK 
APPROACH 


Most TDHF calculations employ Skyrme effective in- 
teraction [57]. The Skyrme energy functional provides 
a good description of nuclear ground state properties 
and collective excitations. Various parametrizations of 
Skyrme effective interaction have been fitted with a dif- 
ferent emphasis on nuclear structure properties. There 
are no free parameters adjusted on the reaction dynam- 
ics. In order to see the parameter dependence, three 
Skyrme parametrizations of SLy4 [58], SkM* [59], and 
UNEDF1 [60] have been implemented to study the dissi- 
pation dynamics in heavy-ion collisions in this work. We 
will not show the calculations with other parametriza- 
tions since they do not produce a significant difference 
for the purpose of the present work. 

The time evolution of self-consistent mean-field gener- 
ated by all the particles is described by the TDHF equa- 
tion 


= 


ihOe = hlp, T, Ë, J, T Ya, (1) 


where the time-dependent single-particle Hamiltonian h 
consists of density p, kinetic density 7, spin density ¢, 
current J, and spin-orbit density J [61]. The main ap- 
proximation of TDHF theory is to treat the many-body 
wave function as an independent particle state at any 
time. Starting from a proper initial state obtained from 
the static HF equation, the TDHF equation is solved to 
determine the wave functions at each time of dynamical 
evolution. The mean-field Hamiltonian h is derived from 
the energy density functional (EDF) 


-> > 


E= I TAT AA (2) 


by using the time-dependent variation 


ðE 
Opji 


hij = (3) 


The energy functional H consists of free kinetic energy, 
Skyrme interaction, and Coulomb energy with exchange 
in the Slater approximation. Note that the TDHF ap- 
proach treats the static properties and reaction dynamics 
in a unified theoretical framework and the same energy 
functional. 

For the sake of later discussions, the time-even and 
time-odd terms of spin-orbit force 


even 1 T T 
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are expressed in terms of various densities and Skyrme 
parameter t4. The index q denotes protons and neu- 
trons. The inclusion of both time-even and time-odd 
terms guarantees the Galilean invariance which should be 
met in a meaningful theory of heavy-ion collisions [61]. 
The odd-odd terms of spin-orbit force are important to 
assure the energy conservation in the free translation of 
a nucleus over the grids as reported in Ref. [42]. 

The set of nonlinear TDHF equations is solved on 
a three-dimensional Cartesian coordinate-space without 
any symmetry restrictions. The derivative is solved with 
fast Fourier transformation in Fourier representation. 
The numerical box of coordinate space is 32x 24x24 fm 
for the present calculations of the light system 1860+160. 
The colliding nuclei locate at an initial center of mass 
distance of 16 fm. Before reaching this distance, the col- 
liding nuclei follow the Rutherford trajectory. The grid 
spacing is taken as 1 fm and a time step of dynamical 
evolution as 0.2 fm/c. The dynamical unitary propaga- 
tor is expanded up to the sixth order of Taylor expansion. 
The choice of these parameters guarantees good numer- 
ical accuracy during the dynamical evolution for all the 
cases studied here. The shift of particle number and to- 
tal energy is less than, respectively, 0.01 and 0.1 MeV 
during the time evolution. 

Fusion occurs when the collective kinetic energy is en- 
tirely converted into the internal excitation of a well- 
defined compound nucleus. TDHF fusion cross section 
is calculated at each energy by the sharp-cutoff approxi- 
mation [62] 


wh? 


D CD min 1 ; $ 
2B ean. deans SY 


Ofus = nis F Ug = 
where u is the reduced mass of the system and Ee.m. is 
the initial center-of-mass energy. The quantities lmax and 
lmin denote the maximum and minimum orbital angular 
momentum for which fusion happens. At low collision 
energy, the minimum orbital angular momentum is zero. 
As the increase of incident energy, there usually appears 
the nonzero lower limit of orbital angular momentum due 
to the transparency behavior in central collisions [36, 63, 
64]. 


II. RESULTS AND DISCUSSIONS 


We have performed fully three-dimensional TDHF cal- 
culations using Sky3D code [65] for heavy-ion collisions 
of 1©O+!6O with modern Skyrme EDF and without any 
symmetry restrictions. This system has been studied 
with TDHF by several groups on fusion, deep-inelastic 
collision, the associated dissipation mechanisms, and the 
extensive comparison with experimental data [8, 9, 12, 
13, 15, 18, 29, 32, 36-38, 42-44, 46, 62]. The present 
studies will be compared with the earlier calculations and 
the available experimental data in order to clarify the 
impact of new terms in Skyrme EDF and the geometric 
symmetries on the reaction dynamics. 
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FIG. 1. (Color online) Final relative kinetic energy as a 


function of center-of-mass energy for head-on collisions of 
1604+160. The solid lines show the results of present studies 
using three Skyrme parametrizations, and the dashed lines 
are taken from Ref. [38] for comparison. 


One of the most interesting experimental observables 
in heavy-ion scattering is the relative kinetic energy of 
the separating ions. Figure 1 shows the final relative ki- 
netic energy as a function of initial center-of-mass (c.m. ) 
energy for head-on collisions of *O+1°O. The solid lines 
with different symbols denote the present studies using 
three Skyrme parametrizations SLy4, SkM*, and UN- 
EDF1. For comparison the dashed lines from earlier 
TDHF calculations [38] are also displayed. The difference 
between present studies and earlier calculations mainly 
lies in three points. First, the axial symmetry was as- 
sumed in earlier calculations, while the present study 
does not impose any symmetry restrictions. Second, we 
include the full spin-orbit (1*s) force, and yet only time- 
even terms of l*s force were incorporated in earlier stud- 
ies. Third, the terms involving the gradient of density in 
Skyrme EDF were replaced with the finite-range Yukawa- 
folding form in earlier studies to simplify the numerical 
calculations, compared with the full treatment in present 
studies. 


The heavy-ion scattering behavior studied here is for 
energies from the threshold of inelastic scattering to an 
initial c.m. energy of 200 MeV. It should be noted that 
TDHF dynamics can only be used to discuss the one- 
body dissipative mechanism. In the energy range of the 
present study, the two-body collision terms are likely to 
play a significant role. However, as presented in the In- 
troduction, beyond TDHF calculations still remain a nu- 
merical difficulty. We then present this work as a study 
of one-body dissipation dynamics. 

As seen in Fig. 1, the overall trend of solid lines is sim- 
ilar for the three parametrizations, and yet there gradu- 
ally appears a distinct discrepancy as the increase of c.m. 
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FIG. 2. (Color online) Percentage of energy dissipation as 
a function of center-of-mass energy for head-on collisions of 
1604+160 with parametrization SLy4. The black, red, and 
blue lines represent the TDHF calculations involving full 1*s, 
time-even l*s, and no l*s force. 


energy. In the energy range studied here, the outgoing 
ions carry the largest kinetic energies using SLy4, and the 
smallest with the newly fitted parametrization UNEDF1. 
We may therefore expect that the energy dissipation will 
be the most distinguished with UNEDF1, and the weak- 
est using SLy4 among the three parametrizations. The 
dashed lines from the earlier studies are shifted notably 
depending on the Skyrme parametrizations. By compar- 
ing the results of the present and earlier studies, the sensi- 
tivity of earlier calculations on Skyrme parametrizations 
is attributed to the approximations stated above on the 
effective interaction and geometric symmetry. It should 
be noted that the final energy has an uncertainty of 2-3 
MeV because the nucleons emitted from the highly ex- 
cited fragments are reflected from the boundaries of the 
numerical box, as seen in deep-inelastic collisions [44] and 
giant resonance calculations [66]. This will not cause an 
essential distinction for the purpose of present studies. 
In the following, the results with SLy4 will be taken as 
an example for further discussions. 


In order to understand the systematic variation of en- 
ergy dissipation as the increase of incident energy, we de- 
fine a quantity which measures the extent of energy dissi- 
pation in deep-inelastic collisions as Pais = 1— Fgn/Ecm., 
where Efn and Ee.m. denote, respectively, the final and 
initial kinetic energy already shown in Fig. 1. The per- 
centage of energy dissipation as a function of incident en- 
ergy is plotted in Fig. 2 for head-on collisions of 1680+160 
with parametrization SLy4. The black, red, and blue 
lines represent the TDHF calculations involving full 1*s, 
time-even l*s, and no l*s force in Skyrme EDF. As can 
be seen, the inclusion of l*s force results in a significant 
enhancement of energy dissipation in deep-inelastic colli- 


FIG. 3. (Color online) Density profile of the separating ions 
at the same relative distance R = 8.3 fm for three incident 
energies Ec.m. = 90 MeV (top), 130 MeV (middle), and 170 
MeV (bottom) for the head-on collisions of ‘*O+1°O with the 
full 1*s force. 


sions. For all three cases the percentage of energy dissipa- 
tion decreases as the incident energy increases. To clar- 
ify the dissipation mechanism, Fig. 3 displays the density 
distribution of separating ions at the same relative dis- 
tance of 8.3 fm for three incident energies of 90, 130, and 
170 MeV in the upper, middle, and bottom panels. The 
three energies locate around the start, middle, and end 
point of the black line shown in Fig. 2. As the increase 
of incident energy the nucleon density around the neck 
region becomes lower and at c.m. energy of 170 MeV the 
nuclei are essentially spherical without a neck formation, 
indicating that fewer nucleons have been transferred be- 
tween the colliding partners. The dynamical mechanism 
in the exit channel of heavy-ion scattering visualized in 
Fig. 3 may be interpreted by the concept of dissipative 
diabatic dynamics which was introduced in the 1980s [67] 
and later applied in the fusion of heavy nuclei [68]. As 
the increase of incident energy, the collective motion is 
so fast that the mean-field does not have enough time 
to rearrange itself and has to keep its identity as much 
as possible. This energy dependence of density evolution 
has also been presented in the entrance channel of fu- 
sion reactions at lower energies for the same and heavier 
systems [13, 14, 50]. 

The percentage of energy dissipation arising solely 
from the time-even and time-odd contributions of l*s 
force are shown in Fig. 4. The energy dissipation caused 
by the time-odd l*s force is equal to the percentage of 
energy dissipation with TDHF calculations involving full 
l*s force minus that with time-even l*s, which can be 
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FIG. 4. (Color online) Percentage of energy dissipation caused 
by the time-odd and time-even contributions of l*s force as 
a function of center-of-mass energy for head-on collisions of 
1604+160 with parametrization SLy4. 


extracted from the black and red lines in Fig. 2. The dis- 
sipation from the time-even contribution is the difference 
between P4is with even l*s and without l*s. As seen in 
Fig. 4, the time-even l*s force provides more important 
contributions to the dissipation than the time-odd terms 
when the relative c.m. energy is less than 160 MeV for 
the '©O+1°O collisions. At low energies the time-even l*s 
force plays a dominant role and the effect of time-odd 1*s 
becomes more pronounced at high energies. This might 
be the reason that earlier TDHF calculations at low col- 
lision energy [37], though the time-odd 1*s force was ne- 
glected, can account for, to some extent, the experimen- 
tal data [69]. The dissipation from the time-odd contri- 
bution increases as the incident energy, as we expected, 
because the current density appearing in the time-odd 
l*s functional is proportional to the velocity for the case 
of a nucleus moving with constant velocity. 


A quantity which expresses the proportion of dissi- 

ae: energy coming from l*s force is defined as Ps.o, = 
(no lxs) (full lxs) (no ls) (full MS) a 

Lee aie where Pj. and Pi 
respond to the blue and black lines in Fig, 2. The aa 
are displayed in Fig. 5 for head-on collisions of *O+1°O 
with the parametrizations SkM*, SLy4, and UNEDF1. 
The percentage of dissipated energy coming from 1*s force 
has similar trend for SkM* and SLy4, but evident distinc- 
tion appears for UNEDF1. As the increase of incident 
energy, the dissipation from |*s force increases and then 
start to decrease after reaching the maximum value of 
65% for SkM*, 63% for SLy4, and 48% for UNEDF1. 
The UNEDF1 calculations predict an overall smaller dis- 
sipation arising from |*s force and lower incident energy 
at the peak of dissipation due to the largest total dissi- 
pation (denominator in the definition of Ps.o.) as seen in 
Fig. 1 and the smallest dissipated energy from |*s force 


(numerator in P;..). The new UNEDF1 parametrization 
was obtained by additional fitting to the experimental ex- 
citation energies of fission isomers in the actinides and is 
suitable for studies of strongly deformed nuclei [60]. The 
origin of the discrepancy between UNEDF1 results and 
those with SkM* and SLy4 parameter sets is still un- 
clear. In spite of the sensitivity of force parameters, it is 
impressed that more than half of the total dissipation for 
SkM* and SLy4, and slightly less than half for UNEDF1 
are attributed to the l*s force. 


We also performed fusion calculations for the 1600+160 
system at a c.m. energy of 70.5 MeV. The reason for 
choosing this collision energy is due to the availability of 
experimental data at the energy range studied here. The 
fusion cross section using three Skyrme parametrizations 
at a c.m. energy of 70.5 MeV is listed in Table I, to- 
gether the available experimental cross section with er- 
rors [70]. TDHF calculations with SkM* and UNEDF1 
parametrizations predict the central fusion collisions at 
the energy of 70.5 MeV because the high-energy fusion 
threshold, as shown in Fig. 1, is higher than 70.5 MeV. 
For SLy4 calculations the central transparency predicted 
by early TDHF calculations [36, 63, 64] appears and the 
lower limit of angular momentum will become nonzero. 
Hence, in the SLy4 calculation of fusion cross section, we 
search both the maximum and minimum impact param- 
eter bmax and bmin at which fusion happens. Here the 
fusion is defined as the compound system evolves long 
enough to offer convincing evidence that the system will 
not undergo separation. The fusion window for angular 
momentum is found to be bmin=1.2 fm and bmax=6.5 fm 
for SLy4, bmax=6.5 fm for SkM*, and 6.4 fm for UN- 
EDF1 within a precision of 0.1 fm. The fusion cross 
section obtained with the three parametrizations over- 
estimates the experimental value by about 20%, which is 
a reasonable agreement considering that no free param- 
eters are adjusted for the reaction dynamics in TDHF 
calculations. The cross section for transparency at 70.5 
MeV is predicted to be quite small with 45 mb for SLy4. 
Whether this effect is real or not would not change much 
the conclusion that the fusion cross section is reasonably 
reproduced by TDHF calculations at this energy. 


The heavy-ion fusion at energies of several times the 
barrier height is a challenge both for theoretical and ex- 
perimental studies. As shown by Simenel et al. [12], for 
the 160+!60 system at Ee.m.=35 MeV, the compound 
nucleus cannot achieve a complete thermal equilibration 
between collective motion and internal excitation after 
the long-time evolution and finally settles at an energy 
difference of about 10 MeV between the c.m. energy and 
the excitation energy. Due to the large remaining col- 
lective kinetic energy, the compound system is likely to 
undergo various breakup and fission events which are be- 
yond the description of TDHF theory. Since these events 
are regarded as fusion in TDHF, a larger fusion cross sec- 
tion is expected at high energies in TDHF calculations. 
In addition, two-body collisions which are likely to play 
an important role at high energies have been neglected 
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FIG. 5. (Color online) Percentage of dissipated energy arising 
from l*s force as a function of center-of-mass energy for head- 
on collisions of *°O+!°O with parametrizations SkM*, SLy4, 
and UNEDF1. 


TABLE I. Calculations of fusion cross section for *°O+!°O at 
Em. = 70.5 MeV with three Skyrme parametrizations and 
experimental data with errors [70]. 


Force fus (mb) 
SLy4 1282 
SkM* 1287 
UNEDF1 1327 
Expt. 1056 + 125 


in mean-field dynamics. From the experimental point of 
view, the fusion cross sections for the 1604+160 system 
have been measured by several groups using different ex- 
perimental approaches and techniques [71-73]. The dif- 
ferent sets of experimental data are not consistent with 
each other at energies of two to four times barrier height, 
leading to a large uncertainty on the experimental cross 
sections [12]. The experimental cross section at c.m. en- 
ergy of 70.5 MeV, although only one set data is available, 
is expected to have large variations as the cases at lower 
energies. Reliable experimental data at well above the 
barrier energy are highly desired. 

The dissipation dynamics has also been investigated for 
medium-mass system of 4°Ca+*°Ca in order to see the 
systematic variations. There is no qualitative distinction 
compared with the light system of '*O+!'6O. Therefore, 
the results for heavier system will not be incorporated in 
this article. 
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IV. SUMMARY 


We have investigated the dissipation dynamics in 
heavy-ion collisions of '©O+!%O using the modern TDHF 
approach. The numerical calculations have been im- 
plemented in a fully three-dimensional coordinate space 
with the modern Skyrme energy functional and with- 
out any symmetry restrictions. The dissipation dynam- 
ics exhibits a universal behavior by using three Skyrme 
parametrizations SkM*, SLy4, and UNEDF1. We re- 
vealed that the percentage of energy dissipation in deep- 
inelastic collisions of the light systems decreases as the 
increase of incident energy due to the interplay between 
the collective motion and the single particle degrees 
of freedom. TDHF calculations with the newly fitted 
parametrization UNEDF1 predict the largest energy dis- 
sipation, while the smallest with SLy4 among the three 
parametrizations. The energy dissipation in the present 
studies has been compared with the results in earlier cal- 
culations. Special attention is paid to the role of time- 
even and time-odd spin-orbit force. The spin-orbit force 
significantly enhances the energy dissipation. The time- 
even coupling of spin-orbit force plays a dominant role at 
low energies, while the effect of time-odd terms becomes 
more pronounced at high energies. More than half of the 
total dissipation for SkM* and SLy4, and slightly less 
than half for UNEDF1 are predicted to come from spin- 
orbit force. We also performed the fusion calculations 
for this system. The theoretical fusion cross section ob- 
tained from the parameter-free TDHF approach agrees 
reasonably well with the experimental data. 
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